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^ ■ Abstract 

. Cumulants represent a natural language for expressing macroscopic proper- 

ties of a solid. We show that cumulants are subject to a nontrivial geometry. 
^ . This geometry provides an intuitive understanding of a number of cumulant 

. relations which had been obtained so far by using algebraic considerations. 

. We give general expressions for their infinitesimal and finite transformations 

0^ ■ and represent a cumulant wave operator through an integration over a path 



X 
J3 



in the Hilbert space. Cases are investigated where this integration can be 



S 

. done exactly. An expression of the ground-state wavefunction in terms of the 

O ■ cumulant wave operator is derived. In the second part of the article we derive 



the cumulant counterpart of Faddeev's equations and show its connection to 
the method of increments. 

I. INTRODUCTION 

In statistical physics, quantum mechanics, quantum field theory and many other fields 
of modern theoretical physics we constantly face objects of the type 

F = ln(e^i^i+^^^^+-). (1) 

Here Ai, y42, ... are abstract elements, Ai, A2, ... are numbers, (...) means some kind of 
averaging and e*- ' -* denotes a generalized exponential function. 
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Assuming the average of 1 to be nonzero, we conclude that F is an analytic function 
of the Aj in a vicinity of zero. Therefore, it has a Taylor expansion with respect to the 
Aj. This expansion is called cumulant expansion, the coefficients of the expansion are called 
cumulants. 

Cumulants (semiinvariants) have been known long before in mathematical statistics and 
probability theory. Kubo ^ pioneered in demonstrating their usefulness in diverse branches 
of physics. Considered as abstract objects cumulants form a special kind of noncommuting 
algebra. In his original paper Kubo applied relations found in this algebra to a diverse 
set of physical problems. In particular, he gave an elegant derivation of the Ursell-Mayer 
expansion for classical and quantum gases that is usually obtained by much longer diagram 
considerations, he rederived Goldstone perturbation theory for the ground-state of a Hamil- 
ton operator, and investigated the motion of an oscillator with a random frequency. Kubo 
stressed the idea that "the generalized cumulant expansion provides us with a point of view 
from which many existent methods in quantum mechanics and statistical physics can be 
unified" [|. 

Cumulants found a fundamental application in solid-state physics (and references 
therein) , [0] , [Q , |^ , as a universal language for expressing properties of a solid in a size- 
consistent way. Expressions for the ground-state energy, expectation values of operators, 
correlation functions and other quantities characterizing a solid were derived by applying 
the cumulant technique. These expressions remain valid in the strongly correlated case, 
when the ground-state of the unperturbed Hamiltonian is not a single particle one and 
the usual diagrammatic theory based on Wick theorem either breaks down or becomes 
very complicated. Also, an important feature is that any cumulant expression vanishes for 
statistically independent processes. Thus the cumulant technique takes into account only 
relevant, statistically connected terms, e.g., of perturbation theory. 

In this article we continue exploring the cumulant technique and show that cumulants 
possess a nontrivial geometry. This geometry helps to make clear and intuitive results 
obtained before by using algebraic considerations. Geometry in an abstract sense means 
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studying how the properties of objects change under transformations 0. First, we find a 
formula for an infinitesimal transformation of a cumulant, and then integrate it in order to 
obtain a finite transformation. 

It comes out that a transformation of a cumulant expression corresponds to a multi- 
plication by some cumulant wave operator Q. Different analytical expressions for Q were 
obtained before through cumulant algebra P| . We give a general expression for Q through an 
integration over a path in Hilbert space, and show cases when this integration can be done 
analytically. Then we investigate the ground-state wavefunction of a Hamilton operator and 
give an expression for it by applying cumulant technique. For this we first need to add a new 
element to the cumulant algebra which is the cumulant of the number 1. Next we show how 
the cumulant expressions found before can be rederived by using geometrical considerations. 
This includes orthogonality relations for Q, for the ground-state energy, and for expectation 
values of operators and correlation functions. Then we give some new applications of the 
cumulant technique. 

First we show how improved short-time propagators for quantum Monte Carlo meth- 
ods can be derived by using cumulant expansions and derive high-order approximations 
for coordinate-momentum and coordinate-coordinate propagators. Then we consider the 
foundation of the method of increments and extend this method to excited states of a 
Hamilton operator. The method of increments has been proven a powerful and intuitive 
way of taking into account electron correlations in solids 0, 0, [|1^, [|ll|], [|1^]. The 



existing theoretical justification of the method of increments is based on Bethe-Goldstone 
equations. The n-th order Bethe-Goldstone equation is equivalent to a variational calcu- 
lation with some trial wavefunction (for the explicit construction of this function see [0). 
This makes the Bethe-Goldstone equations approximate by nature. Therefore, although the 
Bethe-Goldstone equations provide a justification for the method of increments, they are not 
able to predict what part of the problem is neglected by performing an n-th order increment 
calculation. In other words, we do not have an exact n-th order equation, such that the n-th 
order increment equation would be an approximation to this equation. Also we should make 



a remark that the n-th order Bethe-Goldstone equation is restricted to a vacuum taken as 
the unperturbed state. Therefore, it is not apphcable when we need to take into account 
additional correlations starting from an already correlated state. 

On the other hand Faddeev's equations for the 3-body problem are well known |T6 
and have been successfully applied in nuclear physics []T^ and in solid-state theory 



HH) ^3) ^H- Faddeev's equations are exact, this feature makes their application mathe- 



matically controllable. We show that we are able to generate a hierarchy of Faddeev-like 
cumulant equations. Taking the equation of the n-th order in the leading approximation 
we exactly recover the n-th order increment expansion. We can see explicitly from the n-th 
order cumulant Faddeev's equation what part is neglected by the increment approximation. 
Therefore, we generally predict the validity and the precision of the n-th order increment 
calculation. The cumulant Faddeev's equations enable us to formulate an increment method 
for excited states of a Hamilton operator. 

Finally, we give the conclusions and point out future directions of research. 



II. GEOMETRICAL PROPERTIES OF CUMULANTS 

A. Definitions and properties 

We start by defining the cumulant of operators [Q. 
Definition 1 The cumulant of N operators is defined as 

{A,...Ar,r = ^...^ In {x I e'^^\..e'-^- \ y). (2) 

Here x, y are vectors belonging to the Hilbert space, we suppose {x \ y) 0. 

We denote this cumulant as {Ai...A]^y or, when we need to give the averaging vectors 
X and y explicitly, as {x \ Ai...An \ yY . 

Cumulants posses a number of properties (see 0). We briefly repeat them below. To 
simplify notations we consider the case of a two-operator cumulant. 
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1. Linearity: 

{A{aB + PC)Y = a{ABy + P{ACY. 

2. Independence of the norm of the averaging vectors: 

{ex I AB I yy ^ {x \ AB \ y)" for e ^ 0. 

3. Operators A and B are called statistically independent with respect to x,y when the 
average of their product is equal to the product of their averages, i.e., 

{x I AB \ y) = {x \ A \ y) {x \ B \ y) , {x \ y) — 1. 

Then 

{x\AB\ yY = 0. 

4. Applying Definition 1 we can obtain explicit formulas for cumulants. For the case 
{x \y) = 1 they read as 

{x\A\yY = {x\A\y) 

{x \ AB \ yf ^ {x \ AB \ y) - {x \ A \ y){x \ B \ y) , etc. 

5. We need to distinguish the cumulant of two operators A and B and the cumulant of 

an operator AB considered as an entity, when cumulants are evaluated. The latter one we 
denote as {{AB)*). 

Also we need to distinguish between the number 1 and the unit operator 1. For instance 
{l-AY = {AY and {1AY = 0. 

We define formal powers as {A^Y — ( A ■ y Then the formal algebra of power series 

n 

must be considered. 

Let us formulate an important statement. 
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Statement 1 Let H be a Hamilton operator. Let ip be an eigenvector of this operator with 
an eigenvalue E. Then for any vector i/jq, {i/jq \ ip) ^ and for any cumulant series A, not 
including c-numbers the following relations hold 

{iJo\H\ = E (3) 
(^0 I AH I = 0. (4) 

Proof: 

The proof of @ comes from the definition of a cumulant of one operator. Equation 
(I) follows from the fact, that H acting on | ip) leads to a multiplication by a number and 
therefore the logarithm in factorizes. 

The opposite statement is also true: 

Statement 2 Let 

{^0 \AH\^y = o 

for any tpQ, {tpo \ ip) and for any A not including c-numbers. 
Then ip is an eigenvector of H. 
Proof: 

Let P be a projecting operator on the direction of ip. Without lost of generality we 
assume {ipo | ■?/') = 1. Let us consider the following cumulant: 

(^0 I PH I = (^0 \ PH \ ^) - {^0 \ P \ ^){i^o \ H \ ij) = {iJo \ PH - H \ 

Suppose ip is not an eigenstate of H. Then {PH — H) ^ {]. Therefore we can choose a 
vector ipQ such that {ipQ \ PH \ ipY ^ 0. But this contradicts the condition of the theorem. 
Therefore, ip is an eigenstate of H. 
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B. Transformations of averaging vectors 



Here we analyse the basic geometric properties of cumulants. We show that a transfor- 
mation of averaging vectors corresponds to a multiphcation by the cumulant wave operator 
Q (for a definition of Q in terms of power series see [^). 

First we consider an infinitesimal transformation of one of the averaging vectors. 

Statement 3 Let 

= e^^ipB' = (1 + SS)ipB', SS = eS, e « 1. 

Then 

{^a\X\ tfjBy = {^A I X{1 + 5S) I ^B')'. (5) 

Here X is any formal cumulant series. 
Proof: 

The proof follows from the definition of a cumulant. We present it for the case that X 
is a one-operator cumulant. 

(^A I X I ^b)' = |r In {^A I e^^ | e^^V^^') = 
c'Aa=o 

= — In {^A I e^^ I ^B') +e— ttt In {^a \ e^^e^^ | ^b') = 
OA A=o oe e=0 OA A=0 

= {^A\x{i + ss)\^B'y. 

Here we considered an infinitesimal transformation of the left averaging vector. In an 
analogous way, if for the right averaging vector 

= e-'^^A' = (1 + SS)iIja', SS = eS, e « 1 

then (V^ I X I ^b)' = i^ljA' I (1 + 6S^)X \ ^b)"- 
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We conclude that an infinitesimal transformation of the averaging vector result in a 
multiplication of the cumulant operator by 1 + 6S. In the following we call this operator 
infinitesimal cumulant wave operator. 

The next step is going over to a finite transformation by integrating an infinitesimal one. 

We consider the following construction: 

Construction 1 Let ipA, 'ipB, V'c vectors of Hilbert space. Furthermore, let us consider 
a continuous path in Hilbert space between vectors ipB (^nd ipc- Let us choose + 1 points 
ipO) V'l; ■■■ '4'n on this path such that ipQ = ipc o-nd ip^- = ipB o-nd N operators 6 Si such that 

By assuming that the operators SSi are small and using (|^) we obtain 

TV 
i=l 

Taking formally the limit — oo we obtain 

{^a\x\ ^jbY = {i^A I xn I 

N 

n= lim n(l + '^'5i). (6) 

1=1 

Thus we conclude that a transformation of an averaging vector in a cumulant expression 
corresponds to a multiplication with the operator Q. In the following we call this operator 
cumulant wave operator. 

We note that there is a considerable freedom in the Construction 1 for choosing the 
path between vectors ipB, i'c and therefore the operators 6 Si. Thus the operator Q is not 
defined uniquely. Two different representations for Q differ by a cumulant operator with a 
zero contribution to any cumulant expression (such operators can be obtained by integrating 
over a closed path in Hilbert space and subtracting the number 1). 

We also note that Q always contains the number 1 as the first term of the formal series. 

The formula for the finite transformation of the left averaging vector is obtained in 
analogy to (P). In that case we have 
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(7) 



C. The cumulant of the number 1 

To complete the definition of the algebra of cumulant series we need to define the cu- 
mulant of the number 1. By formally putting the number of lambdas in the Definition 1 to 
zero we have 

{x\l\yy = \n{x I y). (8) 

The usefulness of this definition comes from the point that the new object we just intro- 
duced satisfies the general expression for an infinitesimal transformation of a cumulant. It 
can be directly checked that 

{iiJA\i\iJBy = {iJA\i{i + ss)\^B'y. 

Therefore when doing transformations we need not to discriminate between the cumulant 
of the number 1 and "ordinary" cumulants. 



D. Superoperators 

In parts of this article we use cumulants of superoperators, so we will repeat some 
considerations from 0. A superoperator is an operator acting in the space of all operators, 
i.e. Liouville space. The cumulant of a superoperator D is defined in a way completely 
analogous to the ordinary cumulant. For instance: 

{x\D\yr=^ ln(a;|e^^l|y). (9) 

Here 1 is the unit operator of the Liouville space. 

We see that the only difference is that we need to include the unit operator when defining 
cumulants. 
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For some operator H we introduce the superoperator L: L = [H, ...]. 
Cumulants containing L have the following property ( [p^]) 

Statement 4 For any cumulant expression G and for any operator A it is 

{GLAY = {G{LA)'Y. (10) 

Proof: 

The existence of G does not play any role in the following proof and so we omit G to 
simplify the writing. 

92 



{LAY 



In {x I e^i^e^^A | 

Ai=A2=0 



<9AiA2 

ln(x| (l + AiL)(l + A2A) 

OAi CA2 Ai=A2=0 

In (a; I (1 + AiAaLA) I = — In (x | (1 + ALA) | y) (11) 



9Ai 9A2Ai=A2=0 9Aa=o 

Here we wrote A = A1A2. For the second cumulant we have 

{{LA) Y = ^ ln(x|e^^^|y) = ^ In (x | (1 + ALA) | y). (12) 

OA A=0 OA A=0 

So the left and right sides of ( PUD are the same. 

The superoperator L has the following property: if the right averaging vector in a cumu- 
lant is an eigenvector of H than L is equivalent to H in the cumulant expression. Indeed, 
L implies a commutation with H, and AH gives zero, because the right averaging vector is 
an eigenvector of H. 



E. Finite Temperatures 

The above results can be extended to finite temperatures. Here the averaging implies 
taking a trace with some other operator (averaging operator). Therefore, we define the 
cumulant of a product of operators as: 
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{A^...A^ I Gf = ^-j^ InTr (e'^^\..e'-^-G). (13) 

dXi <9A7VAi=...A^=0 ^ ^ 

Here G is an averaging operator. Considering infinitesimal transformations of G of the type 
G = e'^G' we arrive at formulas which are identical to those of the zero temperature case. 

Let us derive the cumulant expression for the free energy of the system. 

Let he H = Hq + Hi where is an unperturbed Hamiltonian, and Hi is a perturbation. 

Then 

F = -/?lnTr(e-'^^) = -/?(1 | e"'^^)' 

In the first step we have used the definition of a cumulant of the number 1 and then trans- 
formed the averaging operator from ^-/^iL+Hi)^-/3Ho e~^'^°. 
Then 5F = F-Fo = -f3{e-^^^+^'^ - 1 \ e"^^")'. 

This expression is analogous to the closed loop expansion of the free energy (thermodynamic 
potential) in field theory. 

III. APPLICATIONS OF CUMULANT GEOMETRY 

Here we use the results obtained in the previous section to describe properties of a 
quantum system. Most of these results were obtained previously by using algebraic consid- 
erations [@-|^ . Here we rederive them from a geometrical point of view. 

A. Energy eigenvalues, orthogonality relations and expectation values 

Statement 5 Let H be a Hamiltonian operator. Furthermore, letip he an eigenvector of this 
operator with an eigenvalue E. Let ipo be an arbitrary vector of Hilbert space {ipo I "V^) 7^ 0, 
and Q a cumulant wave operator transforming ip to ipQ. Then for any formal series A, not 
including c-numbers, the following equations hold 

(^0 I Hn I ^o)' = E (15) 
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(^0 I AHVt I ^o)' = 0. (16) 

Proof: 

By taking equations (^) and transforming the right averaging vector from \I/ to \&o 
we immediately obtain (pTSD, (|T6|). 

Using Statement 2 and transforming \I/ into \&o we find that the opposite statement is 
also true, i.e., if VL transforms \E' into \E'o and ([l^) holds for any A then \E' is an eigenvector 
of H. 

Next we consider expectation values of different physical quantities. 

Statement 6 Let he a vector, with (\E' | \E') = 1, let \E'o he another vector and let Q he a 
wave operator transforming \E' into \E'o- Then 

I G I = (*o I ^^GQ I *o)'- (17) 

Proof: 

First we note that (\E' | G | \E') = (\E' | G | \E')'^. Transforming the right and left averaging 
vectors from ^' to \E'o we arrive at (|T^). 

B. Explicit forms of the Q operator 

We consider cases when the integration (^) can be done analytically. 

Statement 7 Let \E'_b = e^^^c where S is an operator of Liouville space. 

Then one of the possihle forms for the wave operator transforming ^ b i'^to '^c is 

Q = e^. (18) 

In the last expression should he understood in terms of a series expansion. 
Proof: 
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We take SSi in (j^) as: 

5Si = S/N. 

Then we have for Q 

N 

n = hm T\(l + S/N) = e^. 

1=1 

We have used a representation of the formal exponential series known from algebra. It is 
possible to prove a more general statement: 

Statement 8 Let ^b = eP'^^^'^c where S is an operator of Liouville space, and g(S) is an 
analytic function with g(0)=0. 

Then one of the possible forms for a wave operator transforming to is 

Q = e^(^). (19) 

In the last expression e^^^^ should be understood in terms of a series expansion. 
Proof: 

We use for 6Si in (||) the expansion: 

= g\iS/N)S/N 

Then we can write 



N 



n= hm n 

1=1 



S , f iS 



It is known from the theory of series expansions (see for instance [^) that this expression 
tends to 

n = exp [g{S)] . 

The question arises under which conditions an analytic function f{S) can be represented 
by e^*-'^''. The theorem for Weierstrass products PSI states that a necessary and sufficient 



condition is that function f{s) does not have zeroes in the complex plane. 
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For instance, the function 1 + 5* has a zero in the complex plane and therefore cannot 
be represented by e^*-"^-*. Therefore, we cannot state that if = 1 + 5*^^ then 1] = 1 + 5. 

Let us now rederive an operator Q transforming the unperturbed ground-state | ipo) 
of the Hamiltonian Ho + V to the true ground-state | ip) derived in 0. We assume the 
nondegenerate case. 

First we make a remark that 

lim e-^* I ^o) = e-^*(^ I V-o) \^). (20) 

t — S-OO 

Here E is the ground-state energy of H. We assume that | ip) and | tpo) have nonzero overlap. 
Because cumulants do not depend on the norm of the averaging vectors the prefactor in front 
of I ip) is not important and we can write 

lim e"^* \^Ijo)=\^P)- (21) 

t — *oo 

Then using (|TB|) we write 

= lim I e-^'). (22) 

Taking the Laplace transform with respect to the variable z and multiplying by z to get 
the constant part of the expression we obtain 

= hm I -^--) = lim I 1 + -^-r^H) = lim I 1 + -^—-V). (23) 

This expression coincides with the one 0. In the last line we took into account the fact 
that I ipo) is an eigenvector of Hq. 



C. Wavefunction 

Here we consider the question how to obtain the wavefunction from the Vt operator. 

Statement 9 Let ^ and \l'o he vectors of Hilbert space, and Q be a cumulant wave operator 
transforming \E' to '^q. Then there exists an expression for \E' through \l/o and Vt that is he 
derived below. 
Proof: 
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We use the fact that 

V^^i^O I 1 I = VvI;^ln(vI/o I M/) = =1 ^no.m). 

Here 1 is number 1 and | "^norm) means that | is normahzed by the condition (\l/o | = 1. 
By transforming the left side from \& to \E'o we have: 

^norm = left (^0 | | ^o)"- (24) 

Here left means that we should take the gradient only with respect to the left \&o- 

D. Green Functions 

In order to treat Green functions we need to generalize the definition of cumulants by 
including cases with the averaging operators inside. 

We explain this extended definition by giving an example: 

{uAwB)^ = TT^^ In {x I ue^'^^we^^^ \ y). (25) 

oXi aA2 Ai=A2=o 

This is the cumulant of operators A and B with averaging vectors | x) and | y) outside 
and averaging operators u and w inside. 

The usefulness of this definition is seen by noticing, that we can write a Matsubara 



Green's function [24] the following way (we assume {ijj \ ip) = 1): 

{^P I Qe^^c+e-^^ I = I Qe^^c+ | (26) 

Therefore, we can write a Matsubara Green function as a 2-operator cumulant. We can 
find formulas for transformations of averaging operators which are similar to the formulas 
for transformations of averaging vectors. 

As a result we obtain the following expression for the Laplace transform of Rij 

R^A^) = (^ I (^tC*)'4 I (27) 

Z — Li 

This result coincides with the one obtained in by resummation of a power series. 
In a similar way expressions for higher-order Green functions can be obtained. 
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E. Short-time propagators 



The Diffusion Quantum Monte Carlo technique uses Feynman's path integration for the 
evaluation of physical properties of a system. An important issue is the quality of a short- 
time propagatorhj means of which we obtain the evolution of the system after a small period 
of time. M. Suzuki |^ first proposed to calculate improved short-time propagators using 
cumulant expansions. Below we explicitly calculate improved coordinate-momentum and 
coordinate-coordinate propagators for a general physical system with a Hamilton operator 
H = /2 + V{q). The real-space propagator is defined as a matrix element of the evolution 
operator between two 5-functions representing the eigenstates of the position operator |^ . 



Because two delta functions situated at different points of the real space have zero overlap, 
we are not able to apply cumulants directly to the real space. The idea is then to extend 
the method used by Faddeev to obtain the first-order propagator. First we calculate 



the phase-space propagator, and then by taking the Fourier transform come back to the real 
space. 



The phase-space propagator is defined as 



F{p, q, t) = {p\ e-'^'l^ I q) = exp [{p \ e"*^*/^ | g) J . (28) 

Consider H = p^ /2+V{q). The above expression allows us to evaluate an improved short- 
time propagator up to any power of t. Let us calculate it to order t^. Direct calculation 
gives: 



-ipqh-^ - i(pV2 + V{q))th-^ - h^2pW\q) - V"{q))t^ 



(29) 



F{p,q,t) = (27r)-^/2gxp 

Going over to the (q',q") propagator we have: 

F{q', q", t) = {q" \ e^^* | q') = (2vr)-^/^ J e'^^" F{p, q')dp . (30) 

By inserting (^) and evaluating the integral over p we obtain the short-time real space 
propagator: 
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F{q', q", t) = exp 
'{27iihT) 



2th 



■ (31) 



In a similar manner the next order short-time propagators can be computed. 



Let us discuss the physical meaning of the expression (|3ll) . We have two additional 



terms in the propagator as compared to the usual leading order expression derived in p8 
The term {l/2)t^V'{q') corresponds to the fact that the velocity of the particle following 
a Feynman path is not constant, changing due to the force V'{q') and therefore leading 
to a correction of the kinetic energy term. The meaning of the term l/4tV'(g') will be 
considered later. 

Let us use the above results to obtain a high-temperature expansion for the partition 
function of a quantum system. 



The partition function is defined as p8 



Z = Trexp{-(3H) = J {q \ exp(-/3iJ) | q)dq. (32) 
For small (3 we can use the expansion ([3l| ) setting t = —ij3. Then we have to order 

=^= I exp [-^V{q) - (l/4)/.2/3V"(g)] dq. (33) 

To the first order in (3 we obtain a classical expression for a partition function, this 
means that for high temperatures the average energy of the system is large and we can 
neglect quantisation effects. 

Let us discuss the meaning of the correction 1/ AU^ (3'^V" {q). Suppose the potential V{q) 
has a very high but very narrow peak. Then the region of the peak will be prohibited in a 
classical theory because V{q) is large. On the other hand, for the quantum case there will 
be some probability for the system to be in the region of the peak due to tunnelling. For a 
high and narrow peak V"{q) will be large and negative. Now looking at the expression ( |33D 
we see that l/Ah^(3'^V"{q) has then a sign opposite to (3V{q) so the potential is effectively 
smeared out. Therefore, we see that 1/ Afi? jS'^V" {q) term describes an effective smearing of 
the potential due to the quantum tunnelling. 
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The expansion for a partition function obtained above is a high temperature expansion 
(in powers of (3). We note that this expansion does not coincide with the quasiclassical 
expansion (in powers of h). Ahhough both types of expansions give the classical formula 
as the leading term, the higher-order terms are different. The expansion in powers of h was 
first obtained by Wigner To the second order in h it can be written as: 

Z = , ^ I exp \-(3V{q) - {l/l2)h^(3^V"{q) + {l/2A)h^ (3''V' [qf] dq. (34) 

The formulas for improved short-time propagators can be useful to improve the speed 
and the precision of numerical Quantum Monte Carlo calculations |2^. They could also 
improve analytical calculations that make use of short-time propagators. 

We remark that the above technique applies only when the short-time propagator has a 
well-defined Fourier transform with respect to the coordinate. This is usually satisfied for 
potentials approaching zero or a constant at infinite distance. This condition includes the 
majority of potentials used for Monte Carlo calculations of quantum gases and liquids. The 
case when the propagator cannot be Fourier transformed with respect to the coordinate is 
subject to future considerations. 



IV. METHOD OF INCREMENTS 
A. Derivation of the Method of Increments starting from Faddeev's Equations 

It happens frequently that we deal with a Hamiltonian consisting of different parts and 
that these parts considered separately can be solved analytically or numerically. 

The first example of such a situation is iheS-hody problem. Three particles interact with 
each other by two-particle potentials. Considering each pair of particles separately we face 
a 2-body problem which can be solved analytically. Then the question is how we can use 
the information given by these solutions to construct an approximation to the solution of 
the whole 3-body problem. 
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The second example are quantum chemistry calculations where we need to find an ap- 
proximate solution for a large molecule or a cluster having available numerical calculations 
for its subparts, i.e., chemical bonds. 

For the 3-body problem L. D. Faddeev invented coupled integral equations that connect 
the unknown 3-body scattering matrix to 2-body scattering matrices that can be found 
analytically IT^, [[l^] . 



In quantum chemistry the method of increments is well known, that expresses the 
ground-state energy in terms of a sum of local increments given by sets of 1, 2, 3 etc. 
chemical bonds. 

Below we show that Faddeev' s equations and the method of increments are closely re- 
lated. The increment method can be obtained as a leading order approximation to cumulant 
Faddeev-like equations that we derive below. 

We start by writing the wave operator in the form p[. 



l + -^V). (35) 



\n) = lim _ 

We also define a scattering operator S, which is small for small V, 

Let us first consider the case when the perturbation consists of three parts. 

Statement 10 Consider a Hamiltonian with three perturbations, Hi, H2 and H^ 

H = Ho + Hi + H2 + H3. (37) 

Let Si, S2, S3 be scattering operators for Hq + Hi, Hq + H2 and Hq + if 3 respectively. 

Then the scattering operator S of the Hamiltonian H can be found from the solution of 
the formal equations 

Ti = Si{l + T2 + T,) 

T2 = S2{l + Ti+n) ^ ^ 

(38) 

T3 = S,{l + Ti + T2) 

S = Ti+T2 + Tg. 
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Proof: 

We write in a power series as 



S=E[ TriHi + H, + Hs)) (39) 



and define 



Ai = -^-ftHi, a, = -^-frH2, A, = —^H^. (40) 
z — Ho z — Ho z — Ho 

Then 

oo 

S=E {A^ + A2 + Asy' = EA^... + EA2... + EAs... = T,+T2 + Ts. (41) 

n=l 

We divided the terms in the power series into three groups, the terms in the first group 
begin with Ai, the ones in the second group with A2, in the third group with ^3. We call 
these three sums Ti, T2, T3. 

We rewrite Ti in the form 

Ti = {A, + Aj + ...) + iAi + Aj + ...)A2 + {A, + Aj + ...)As + ...= ^ ^ 

(42) 

= {Ai + Al + ...){1+T2 + Ts). 

The first bracket defines Si. The expressions for T2, T3 are similar. This results in 
Faddeev's equations (^8|) . 

In the original work of Faddeev 5*1, 5*2, S3 correspond to two-body scattering matri- 
ces coming from three pairs of particles in the 3-body problem, and S corresponds to the 
unknown 3-body scattering matrix. 

We easily generalize our equations to the case when the perturbation has parts, in 
that case we have N equations of the type 

T, = S,(^l + ^T}j. (43) 

Assuming Si to be small, the equations can be solved by iterations. Then 5* is given by (pl). 
In the leading approximation we have 
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T, = S,, S = Y.S,. (44) 

By using the formulas E = {H \ Q) and $7 = 1 + S* where E is the ground-state energy 
of H we can write 

5E = E-Eo = J26E,, (45) 
6Ei = {H\S,). (46) 

Here Eq is the ground-state energy of Hq. The equation (|^) gives the simplest example 
of a increment method. We call the SEi^s first-order increments. 

Let us state the result: we start from Faddeev's equations, solve them in the leading 
order approximation and obtain a first-order increment method. 

Below we show that we can formulate higher-order Faddeev's equations. Taking the 
N-th order Faddeev's equation and solving it in the leading order approximation we obtain 
exactly the N-th order increment method. Therefore Faddeev's equations and increment 
methods are closely related. 

As an example we derive the second-order Faddeev's equations. We again take a Hamil- 
tonian consisting of Hq part and perturbations 

H = H0 + H1 + ... + Hr,. (47) 

But now we additionally assume that we can solve analytically or numerically all the Hamil- 
tonians of the form Hq + Hi + Hj, i.e., that we know the scattering operators Sij for these 
Hamiltonians. 

Then we write again the expression for S, 

00 

s=j:{a, + ... + A^r = 

n=l 

EAi...A2... + E^...v43... + E^i...v43... + ... + E^=i^i + ...E:Li^5^= (48) 

This resummation is more sophisticated and needs explanations. First, we open brackets in 
the sum in order to obtain a power series with respect to Aj's. For each term in this power 
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series we look for the first two letters in it, and regroup terms according to this property. 
For example, Ti 2 will be the sum of all terms in the power series having Ai and A2 as the 
first two letters, so that ,e.g., the term A\A\Ai is included in T1.2. 

We are still left with terms containing only one letter like A\. These terms are not 
contained in any of Tj^'s so we sum them separately, this explains the result of 

Now we consider Ti 2 : 



(49) 



Ti,2 = (A1A2 + A,Al + AIA2 + A.A^A, + ...) [I + E(,,,)^(i,2) + ^^^1.2 

= (51,2 - Si- S2) (l + Xl(ij)^(i,2) + ^*7^1,2 Si 

This result is obtained by considerations similar to the previous ones. The first bracket is 
a sum of all possible terms consisting of A^s and y42's only. The second bracket is a sum 
of all possible terms that do not have Ai and A2 as the first two letters. 5*1,2 denotes the 
scattering operator for the Hamiltonian Hq + Hi + H2. 

The similar equations are obtained for all other Tij. Together with (^) these equations 
give the second-order cumulant Faddeev's equations. 

Taking again the leading approximation for the equations we have: 

S = J:S^ + T.K^J (50) 

where 

Kij = Sij Si Sj. (51) 

Here Kij is a cumulant operator corresponding to the second-order increment. 

Using the formula for the ground state energy E = {H \ Q) we arrive at an expression 
up to the second-order increment expression 

5E = Y.5Ei + Y.5Ei,. (52) 

i ij 

The second order increments SEij are given by 

SE,^ = {H\K,,) = {H\S,,) - 6E, - 6E, (53) 
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By continuing this procedure and taking the leading approximation to the next contribution 
we end up with the third-order cumulant expression for S 

S = Y.S^ + T.K^, + Y.L^,k (54) 

i i,j i,j,k 



■'ijk 



Sijk — Kij — Kik — Kjk — Ki — Kj — Kk- (55) 



Again this expression agrees up to third-order with the energy expression of the increment 
method. The same holds true for higher-order contributions. 

Next we consider the error when terms up to the N-th order increment are taken into 



account. From ( ^Qf ) we see that the first bracket is of second order in the perturbation, 
while the remaining term that we omit in the leading-order approximation is of first order, 
therefore making a leading-order approximation to the equations finally leads to an error 5* 
and E which is of third order. Generally, inclusion of Nth order increments will lead to an 
error which is of the order + 1. 

Actually, the increment expansion may converge much more rapidly and the perturbation 
does not need to be small. 

Consider for instance a large molecule. The term omitted by the leading order approx- 
imation to Faddeev's equations can make a significant change to the result only when all 
its parts add coherently; otherwise its contribution will be close to zero. Consider a term 
corresponding to a set of K chemical bonds. If the spatial extent of this set is larger than 
a a characteristic length we called the correlation length A, then the omitted terms in the 
Faddeev's equation will sum incoherently and can be neglected. Therefore we can terminate 
an expansion of the correlation energy in terms of increments when the spatial extent of a 
set of bonds is of the order of the correlation length. This is in accordance with the well 
known idea, that the minimum size to which we can reduce a system not changing its macro- 



scopic properties is dictated by the correlation length In order to formulate it more 
rigorously we may state that the macroscopic properties of an infinite solid can be described 
by an incremental expansion with increments taken up to the correlation length. Therefore 
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incremental methods will work efficiently for systems having a moderate correlation length, 
like semiconductors [ITT . 



B. Method of increments for excited states 

In the following we formulate a method of increments for excited states. We begin with a 
nondegenerate unperturbed Hamiltonian, so perturbations do not split excitation energies, 
instead they only shift them. 

Consider again as an example a Hamiltonian with three perturbations H = Hq + Hi + 
H2 + H3. 

Consider an excited state Ei of Hq. We want to show that the incremental method for 
an excited state can be formulated in a similar way as done before for the ground-state . 
We start from a formula derived in for the excited states of a Hamiltonian: 

E, = I Hn I (56) 

Here | ipi) is an excited state of Hq having an overlap of more than 1/2 with the true 
excited state we are looking for (for details see 

Let us apply the incremental expression for Q to this formula. 
To first order of the incremental expansion (see (PD) we obtain 

5E^ = SEi^i + 5Ei^2 + 5Ei^^. (57) 

Here SEi i,SEi^2,SEi 3 are increments corresponding to the excited state i, 

SEij = I HS, I tfji)'. (58) 

In the second order of the incremental expansion we find: 

6E, = 6Ei^i + 6Ei^2 + 5^^,3 + 5^.,l,2 + 5^i,2,3 + 5^i,l,3 (59) 

where 5-Ej,i,2 = {ipi \ HSi^2 \ 'ipiY - - 5-E'j,2 etc. 

In an analogous way next order increment expressions can be obtained. 
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Let us consider as an example of an increment calculation the valence band calculations 



done in 31 



In a first step self-consistent field CSCF calculations are done and the energy bands in 
SCF approximation are obtained. Then the SCF states are transformed into Wannier states 
(localized orbitals) and the residual part of the interaction if^es = H — Hscf is expressed 
in terms of matrix elements between the localized bonds. 

The increment method up to forth order was used to obtain the valence band. In this case 
SCF hamiltonian plays the role of Hq and the matrix elements of the residual Hamiltonian 
between pairs of localized bonds play role of perturbations. In the calculations increments 
corresponding to NN (nearest neighbours) and NNN (next-nearest neighbours) were retained 
and the effect of more distant correlations was approximated by introducing a polarizable 
continuum. 

The results have been shown to describe with a good accuracy the reduction of the 



valence bonds of diamond due to correlations 31 



V. CONCLUSIONS 

In this article we considered the question, how geometrical methods can be used to 
formulate the theory of cumulant expansions. Cumulant geometry simplifies the proofs of 
many expressions obtained before by cumulant algebra and makes them transparent. 

Then we studied applications of the cumulant technique to short-time propagators and 
to the method of increments. In both cases the use of cumulants is a natural choice. A 
hierarchy of cumulant Faddeev-like equations was generated and a method of increments for 
excited states was established. 

Future directions of research will include diverse applications of the cumulant technique 
to solid state theory. A particularly interesting question there is how properties of the 
infinite solid can be related to properties of a finite cluster. This question is of outmost 
importance for all numerical calculations. 
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